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ABSTRACT 




Solutions to a nonlinear wave equation were analyzed for 
their stability. The wave equation is a Klein-Gordon equa- 
tion with the mass replaced by the square of the wave func- 
tion. This wave equation has propagating solutions which 
are unbounded or periodic, depending on the sign of the non- 
linear term and the propagation speed which can be sub- or 
super-light velocity. The stability of the periodic sub- 
light velocity solution was investigated by the method of 
characteristic exponents and was found to be indifferent. 
Liapounoff's direct method and Sturrock's analysis of the 
dispersion relation combined with a WKB technique were ap- 
plied to a linearized perturbation on a static solution of 

2 

the field equation. The periodic solution with 6 < 1 is 

stable, while the method of characteristic exponents gives 
indifference. The super-light velocity solutions are un- 
stable. Due to the limitations of the approximations, it 
could not be determined whether the instabilitv is absolute 



or convective. 
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INTRODUCTION 



In recent years there has been increased interest in 
nonlinear wave equations in connection with theories of ele- 
mentary particles. Such equations arise when one attempts 
to describe a particle or a system of particles, usinq the 
theory that the particle itself arises from self-interactions 
within the wave function. Heisenberg's nonlinear spinor 
theory, [1], a quantized version of such a theory, was able 
to reproduce certain features of the spectrum of elementary 
particles . 

Nonlinear differential equations have solutions that 
are fundamentally different from the solutions of linear 
equations . [2] They can disappear, become singular, or arise 
from nothing. A central problem lies in the stability of 
the solutions. 

The purpose of this investigation is to study a possible 
method of stability analysis as applied to the solution of a 
wave equation which may serve as representative of a whole 
class of nonlinear wave problems. 

Three methods of stability analysis will be employed. 

Two of these, Liapounoff's method and the method of charac- 
teristic exponents, [3], are interrelated. LiaDounoff's 
method of stability analysis concerns the stability of the 
equilibrium points. A theorem, also due to Liapounoff, 
states that the results regarding stability of the linear- 
ized theory also hold for the nonlinear case, provided 
Liapounoff's method is applicable (the characteristic 
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equation does not have pure imaginary roots) . The method of 
characteristic exponents concerns the stability of periodic 
solutions . 

The third method is the analysis of the dispersion re- 
lation as described by P. A. Sturrock [4], which was devel- 
oped from stability problems in plasma physics. This method 
is also only applicable in the linearized theory. 

The representative differential equation to be used is 

V 2 $ - e\ 2 <t> 5 - - 2 = 0 , 

c at 2 

which will be derived in Chapter I from the Klein-Gordon 
equation . 

The particular questions being asked are: 

1. What are the static solutions of the equation 
for given symmetries? 

2. Are time dependent solutions to perturbations 
stable, or for what wave lengths are they 
unstable? 

Chapter I of the paper contains the derivation of the 
wave equation to be used, an analysis of the types of solu- 
tions and some of the pertinent solutions. In the second 
chapter selected solutions and perturbations are analyzed 
for stability. Chapter III is a summary of results and con- 
clusions drawn from those results. Two appendixes attached 
to the paper contain expansions or alternatives to solutions 
derived in the paper itself. 
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CHAPTER I 



THE NONLINEAR WAVE EQUATION 



1.1 DERIVATION OF THE NONLINEAR EQUATION 

The wave equation chosen as a starting point for the 
derivation is the Klein-Gordon wave equation for a free 
particle: 



Two reasons may be advanced for this choice. First, 
the Klein-Gordon equation is the basic equation for relativ- 
istic quantum mechanics reflecting relativistic energv con- 
servation. Therefore a nonlinear equation based on this 
ought to contain similar features which could be interpreted 
on the basis of the known. Second, both the electromagnetic 
and Yukawa static potentials can be derived from the relativ- 
istic equation by setting the mass term equal to zero or the 
mass of a pi-meson, respectively. Thus, one can hope to find 
in the static solution of the nonlinear theory the analogues 
to the known static potentials. 

To derive a nonlinear wave equation, an approach is 
taken similar to the one that led Heisenberg to his nonlinear 
spinor equation. [1], Instead of making the mass term a con- 
stant, we choose to write m ~ <j> 2 , and replace the entire term 



whe re e = ± 1 . 

This choice appears to be plausible by recalling that 
in quantum mechanics, the probability of finding a particle 




( 1 ) 



by 




J 
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at the position r is P(r)dr = 4>*4>dr. The particle can be 
thought of as being smeared out over soace and the require- 
ment that 



means that the particle remains intact. In a certain sense 
then, the probability of finding a particle at the position 
r could be construed to indicate its mass density at the posi- 
tion r. The wave equation to be used is therefore 



e = 1 corresponds to a modified Klein-Gordon equation with 
the mass replaced by self-interaction of the field. e = -1 
represents the equation of an electromagnetic potential with 
a source term given by the field self-interaction. 

1.2 BEHAVIOR OF THE STATIC SOLUTIONS 

The static potential function is the solution of the 
equation 



The equation can be solved using spherical coordinates 
and assuming that the solution is spherically symmetrical. 
The equation to be solved is therefore: 




v 




( 2 ) 




(3) 




0 * 
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To proceed, assume the solution is of the form 



where u is an unknown function of r. By substitution one ob- 
tains : 



n d 2 F 

: 5? 



(u* ) 



r n g u" + 2 (n + l)r n - 1 g u' 



+ n(n+l)r n 2 F = eX 2 r sn F 5 



One may choose n and u arbitrarily, and then determine F. 
If u = In r, then the differential equation becomes 



n-2 d 2 F n-2 dF , . . . n-2 dF ^ n-2 „ , 2 5 n „ 5 

r r - 3 — + 2 (n+1 ) r -g— + n(n+l)r F = eX 2 r F . 

du 2 du du 



Choose n so that n - 2 = 5n; n = The differential 

equation then reduces to: 

- 4fF - eX 2 F 5 = 0 . 



d 2 F 



du' 



Let 



z 



dF 

du 



dz _ 1 
dF z 



d 2 F 

du 2 



or 




d 2 F 

du 2 



There fore : 



z = eX 2 F 5 + %F. 

z 2 eX 2 F 6 , F 2 _ u 

2 6 8 * 



(4) 



Note that this equation resembles the energy integral of clas- 
sical mechanics where z is analogous to velocity, the function 
of F to a negative potential function, in which F plays the 
role of the space coordinate, and h to the total energy. One 
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can write the equation in the form: 



z 



+ 




( 5 ) 



One can ask what a plot of z versus F (a coordinate- 
velocity phase diagram in the mechanical analoq) would look 
like for various values of h and e, provided one requires that 
z and F remain real. A qualitative analysis yields the re- 
sults in the following sections. 

1.2.1 THE CASE e = 1 
If e = 1, then 



If h > 0, z has a minimum of /2h when F = 0, but other- 
wise behaves much in the same manner as for h = 0. 

If h < 0, z can take on all values, and for large F var- 




If h = 0, then for large F 




ies as ± 




In this case, however, when z is zero the 



absolute value of F is greater than zero. 



That is 



| F min | = a > 0 . 



A little algebra will yield the one real value of a as 



a = ( A+B ) ^ , 



1 



A = jj- [12b 2 A + Vl44bjx 2 +1] 3 



1 



B = rpj- [12b 2 A - Vl44bjA 2 +l] 3 



and 



b 2 = | 2h | . The results are shown graphically 



in Figure 1. 
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z 




FIGURE 1 



Phase Diagram vs. F) for Static Case, e = 1 

On each curve the curve parameter is u. The arrows indicate 
the direction in which the parameter, u, increases. The ac- 
tual values of the parameter along the curves depend on the 
particular solution. 

1.2.2 THE CASE e = -1 

In this case equation (5) becomes 

z = i^h-Alll + %f 2 . 

If h > 0, a phase plot will yield a closed loon about the 
origin, where values of z and F outside the loop are forbidden. 
The characteristics of the loop are dependent on the size of 
h and A. 
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First, one notices that 



= /2h 



when 



F = 0 , ± 



M 4 

4 A ' 



If z is set to zero, one may solve an equivalent cubic 
for F max. 



3 3v 

v - — 



3b 2 
0 



= 0 . 



4 A 2 A 2 



F max = ± (A+B)^ if 144A 2 bJ > 1 



A = [ 12 A b 2 + Vl44A 2 bJ - 1 ] 3 , 



B = JX [12Xb 0 _ V 144 A 2 b; - 1] 1 , 



and 



b 0 = 2h 




FIGURE 2 
dF 

Phase Diagram vs. F) for Static Case, e = -1, h j> 0 



When 144A 2 bg < 1, 



1 X' 

F max = ± [y cos 0/3] * , cos 0 



12b 2 A 
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There is only one root for F since if 144A 2 b>Q > 1, two 



roots of v are imaginary; if 144A b^ _< 1, two roots of v are 
negative and would yield an imaginary F. These results are 
illustrated in Figure 2. 

If h = 0, the double loop in Figure 2 is obtained. This 
case is only a degenerate case of h > 0. The two nodes at 
F = 0 meet at the origin. In this case 



F max = ± 



IT 2 



% 



When h < 0, equation (5) may be written 



+%F 2 , 



where 



z = ±l/-b?. F 



bj = |2h 



F now has both a maximum and a minimum. The absolute value of 
the minima must be greater than zero. The effect is to sep- 
arate the two lobes for h = 0 in Figure 2 and move them along 
the F axis. From the solution for the roots of the cubic, 

3 3 3b o _ 0 

v ‘ IX 2V + x 2 " 0 ' 

it can be seen that there is no real F if 



144A 2 bJ > 1 



because the roots for v are either negative or imaginary. If 

144A 2 bJ < 1 , 

there will be two solutions for F. These represent the maxima 
and minima for F. 

F max, min = ± (^cos 9/3)^ , - s i n 9/3-cos 9/3)] 3 

where cos 9 = -12Abg . 
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If 144A‘bg = 1, there is one real root for F, 

p = + ( — ) ^ 

“ 4a ; 

This represents a point solution for F. The results are di- 
gramed in Figure 3. 




FIGURE 3 



,dF 



Phase Diagram vs. F) for Static Case e = -1 , h < 0 



All the solutions for <)>q will be periodic except when h = 0, 
or when F recedes to a point when h < 0. 



1.2.3 ANALYTIC SOLUTIONS FOR <|> WHEN h = 0 

The two most interesting cases in the last two sections 
are those for which h = 0. 

When e = 1, h = 0, one obtains the equation 

dF, 2 _ A 2 F 6 ^ 2 

W + * 

The equation may be solved by letting 



F = v 



-* 



Then 



3(g) 2 = 3v 2 + 4 A 2 
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^ = In r/K 




which upon transposing the sinh ^ and using the exponential 
transformation for the sinh, One obtains : 




/3<r 




/3 k 



. X ( r 2 - k 2 )_ 



This solution, (see Figure 4) , represents a potential of 
infinite height at r = k. Suppose this potential represented 
the potential from a particle. Such a potential would indicate 
that the particle consisted of a central core of radius k in 
which there is no real potential. Outside the core, the po- 
tential varies much as a classical 1/r potential. The concent 
of a hard core is not new. Arguments have been advanced that 
such a core might exist in nuclei between particles. [5] 




r 



k 



FIGURE 4 



Static Solution (4 >q vs. r) e = 1, h = 0, r > k 
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It is surprising that one gets a long range potential in- 
stead of the short range Yukawa potential of the Klein-Gordon 
equation . 

Using the solution for 4 >q one can determine the applicable 
regions of the phase diagram in Figure 1. 

Since 



F 

F 




+ 



/3<r 2 

A (r 2 -K 2 ) 



When F = 0, r = °°. That is, the saddle point in Figure 1 
belongs to u = ®, 

When considering the derivatives, one obtains 



dF 1 dF 

du 2 " r dr * 

dF 

When F > 0, ^ < 0. Therefore, the applicable region of 

dF 

the phase diagram is the fourth quadrant. When F < 0, ^ > 0. 
The applicable region in the phase diagram is the second quad- 
rant. 



The first and third quadrants must then belong to a dif- 
ferent solution. In arriving at the solution, the equation, 

3 ( 3u } = 3v ' + 4X2 ' 

was derived and the solution was obtained using the positive 
root for the derivative. If now the negative root is used, 
one will eventually get 

-1 v/3 



sihn 



2 A 



= In K/r . 



The solution for <f>g in this case is 



*0 = 1 






A (K 2 -r 2 ) 
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This solution represents a potential well, the general features 
of which are shown in Figure 5. 




X 



K 



FIGURE 5 



Static Solution (<)>q vs. r) e = 1, h = 0, r<K 

This solution is represented in Figure 1 by the first and 
third quadrants. The saddle point belongs to r = 0 , or u = 

If one solves the static equation for e = -1, h = 0, 
quite a different result occurs. 

Let F = v 2 , Then 



where k is the integration constant. Transposing the inverse 
cosh, substituting the exponential equivalent for the cosh, and 
then solving for v, will yield 




2 



3v 2 - 4 A 2 



which has the solution 



, -1 v/T . . 

cosh —fpr = In r/K 



v 



A (r 2 +K 2 ) 
rK/3 



<J> 0 = (rv) 



k/3 



(7) 



A (r 2 +K 2 ) 
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* 

Equation (7) represents an equation which bears a strik- 
ing resemblance to a static electromagnetic field potential. 
Besides the resemblance to coulomb's law, this solution has 
another point of similarity with electric potential fields. 

The differential equation from which it was developed is equa- 
tion (3) with £ = -1, 

V 2 cj> 0 = -X 2 (^ (8) 



One might compare this with the static potential equa- 
tion with a source term, 




where <p represents potential, p the charge density, and e Q the 
permittivity of free space. This potential function, and the 
wave equation it implies, shall be used for the time dependent 
stability analysis in the following sections. Other solutions 
of the static wave equation may be found in Appendix A. 

Assume that the 4 >q of equation (7) represents the coulomb 
potential from an electron. One might ask what would the con- 
stants k and X be, and what would they represent. 

To determine the constants two relations will be required. 

t 

The first is derived from the fact that in electromagnetism 

V 2 <f> = - , 

8 0 



and 




V 2 (pdv = -2— 



where q is charge. In this case q = e, the electronic charge. 

* Chandrasekhar [6] quotes this solution in the theory 
of stellar structure. 
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But for the nonlinear wave equation 



v2 *0 - - 75 

A 



:/!' 



r 2 + < 2 



and in spherical coordinates 
r* 



— = d0 
e 0 

0 Jo 



r2-n r' 
dijj 

0 



5 5 

1 2 r 2 dr sin 0 

A^ 5 



(r 2 + < 2 ) 2 



e_ 



5 5 

A -4 2 

4tt • 3 k 

A* 



r°° 



r dr 



0 (r 2 + < 2 ) 2 



Let r = k tan0. 

Then the integral becomes 




sin 2 0 cos 0 d0 



e_ 

e, 



4ttk^3^ n 1 A 

r — = 9*10 4i\e 

A* 



and 



k 



1 



(£) ‘ = 3 4 ♦ 10 9 e 



'A 

Conversations with those performing the collidine beam ex- 
periment at Stanford University indicate that coulomb's law 
holds to within 2% at r = 3*10 ^ meters. One mav use this 
data to obtain a second relation between « and A. 

From the ratio for k/A evaluated above, 



*0 = 



9 ° 10 9 e 
(r 2 + < 2 ) ^ 
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Now, let ^ and £ 2 defined so that 



5 1 r ' 



(r 2 +K 2 ) * 



A plot of these two functions is shown below in Figure 6 




Comparison of Coulomb Potential and <J>q For e = -1, h = 0 

At r^ , ^ cannot be more than 2% below That is 

- 98 ^i at r o * 

This, in effect, puts an upper bound on k. 

Let 

?2 = • 98 ? 1 at r Q . 



Then : 



§£ • io 17 



[9*10 ^^+k 2 ] ^ 
.-17 



.-17 



k max = .610 

The range of k is 

0 <_ k _< .6*10 

using the upper limit of < to obtain A, one would get 

A = 5. 
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1.3 BEHAVIOR OF THE SOLUTIONS OF THE TIME DEPENDENT EQUATION 



One can find propagating solutions to the time dependent 
equation in one dimension. Such a solution will give some 
idea of the behavior of the solutions for more complicated 
symmetries. The time dependent wave equation in one dimen- 
sion is: 

ill _ 111 = e x 2 5 . 

9x 2 c 2 at 2 



Let 

x = ct , 

S = x - 6 • T , 




Assume that <t> is of the form 

<P = <Ms ) , 



i.e., we confine our attention to propagating solutions. 
Then 



9 2 <)> _ d 2 (j) 
9x 2 dS 2 

and 

111 = 32 d^_ 

9 t 2 dS 2 

The differential equation is 



d 2 <p _ £ A . 2 <j ) 5 
dS 2 1-6 2 



Let 




Then 



and 




As in section 1.2 the first integral yields an equation 
similar to the energy equation in classical mechanics, 



A qualitative analysis of this equation using the methods 
in section 1.2 is presented in the following section. 

1.3.1 DISCUSSION OF \p vs. <p 

If in the equation for ip , 



curves of (f) vs. will be the same as those in Figure 1 with 
F replaced by ij> and z by ip . Again this case represents solu- 
tions which are unbounded. 

If e = 1, B 2 >1, the second term in the equation for \p 
is negative and h must be greater than 0. The phase diagram 
will look something like Figure 7. 



2 6 

6 ( 1 - 3 2 ) 



h 



or 




e = 1 , and B 2 < 1 




FIGURE 7 

Phase Diagram vs. <(>) for Wave Equation With 

e = 1 , 8 2 > 1 , or e = -1, 6 2 < 1 
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One has the result then, that for e = 1, real sub-light 
velocity solutions exist which are unbounded. The only peri- 
odic solutions are those which exist for super-light veloc- 
ities . 

No other cases need be considered, since if e = -1, one 
merely needs to invert B 2 in the last two cases to obtain the 
same results. That is 



1.3.2 ANALYTIC SOLUTIONS OF THE ONE DIMENSIONAL EQUATION 

Because of the redundancy of cases the solutions shall be 
carried out with B 2 < 1 • 

Let 



B 2 > 1, e = -1 corresponds to B 2 < 1, e = 1 , 
and B 2 < 1, e = -1 corresponds to B 2 > 1, e = 1 . 



2h = 6b 2 ; 6 = ±1 . 




( 10 ) 



where 




Let 




then 




Let 




then 
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where the limits on the integral have been left open since 
they depend on e . 

If e = 6 = 1, 
the solution is 



cos 0 = cn(y,k) , 

where cn is an elliptic function with modulus k. The variables 



are 



cos 0 



1 + (l-/3)v 
1 + (1+/3) v 






2 _ 2 + /3 



Y 



2 a o b o 3 * (s - s o 



) . 



Substituting for cos 0 in the solution, one obtains v as 

1-cn y 



v = 



(/3-l) + (/3+l) cn y 



and 



. 1 % 1 

<b = — v = — 

a 0 a 0 



1-cn y 



( /3-1) + ( /3+1) cn Y. 



% 



Let 6 = -1, £ — 1 . 



Then 



r v 



dv 



3 -,M% 



1 [v(v -1) ] 



2a 0 b 0 (S S Q ) 



Substitute v = — 



w 



Then 






w 



dw 



(l-w 4 )^ 



2a 0 b 0^ s S 0 ) 
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The solution is 



w 



= cn [2/2a 0 b 0 (S-S Q ) ] = cn y, 



with modulus 



k 2 = k 



<l> = 



a 0 /cnr 



Then 



One other solution may be obtained when b 2 , or h = 0 . 



must be greater than zero. In this case the equa- 



l-6' 
tion becomes 



d<J> _ X<p 
dS 



2 l 1% 



= a<f> , 



[ 3 | 1-8 | ] 

where a can be either positive or negative. 

<f> = ± [2a (S Q -S ) ]"* . 



These solutions are all unbounded. It might be of inter- 
est to construct a phase diagram similar to that of Figure 1 
with the limiting values of S indicated. The h = 0 case is 
represented by the last solution obtained. When a is positive, 
and <j> is negative, the third quadrant curve is obtained. The 
saddle point belongs to S = in this case. If a is nega- 
tive, the second and fourth quadrant curves are applicable. 

The saddle point belongs to S = 00 . The results are shown in 
Figure 8. The arrows indicate the direction of increasing S. 

When h ^ 0 , the other curves are obtained, in general, 
using the same methods as in section 1.2.1. If h > 0, — is 
a minimum when <}> = 0 . This occurs when cn y = 1 or S = Sq . 

The derivative approaches infinity when <p approaches infinity 
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FIGURE 8 

Phase Diagram (^| vs. <|>) for Wave Equation with e= +1, 3 2 < 1 



or when 

/3-1 = ( /3+1) cn y , 

for the solution obtained. The cn y lies somewhere in its 
second quarter cycle. That is 

K < y < 2K , 

where K is the complete elliptic integral of the first kind. 

S would have some corresponding value, , as labeled on Fig- 

ure 8 . When h < 0, 33 - = 0 when £ = 1. This corresnonds to 

dS 

cn y = 1, or S = S^. The derivative approaches infinity when 

cn y = 0 . 

This occurs when cn y is at its quarter cycle, or S = S^. . 
This is the label used in Figure 8 . 

These solutions became infinite for certain values of S. 
Adjacent to these infinities the function is not real. It is 
therefore not possible to construct at any time, t, a 
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perturbation which is real for all space. These solutions do 
not allow a ready physical interpretation. 

If 5 = 1, e = -1, the equation becomes 



dv 



[v(l-v 3 ) ] 



77T¥ “ 2a o b o (s - s o } 



The solution is 






cos 0 = cn (2aQbg3 4 (S-Sq) = cn y , 



where the modulus of the elliptic function is 

k 2 = 2^3 

K -3 



and 



cos 0 = izliiHiv 
l+(/5-l)v 



Upon making the proper substitutions and solving for v, 
one obtains <)> ; 




1-cn y 

1+/3+ ( /3-1) cn y 




( 11 ) 



This solution represents the only type of solution that 
is oscillatory in nature. 
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CHAPTER II 



STABILITY OF SOLUTIONS TO THE WAVE EQUATION 



In this section, the stability of the solutions to the 
wave equation is investigated. The time dependent periodic 
solution is investigated using the method of characteristic 
exponents. The following sections will discuss the stability 
of wave-like perturbations on the static solution for e = -1. 
The methods to be used are the characteristic exponent method, 
Liapounoff's method and Sturrock's stability analysis. Each 
of these methods shall be described as they occur. 

2.1 THE UNBOUNDED TIME DEPENDENT SOLUTIONS 

The solutions obtained in section 1.3.2 with e = 1 in 
equation 10 all yield unbounded functions, not real for all 
r and t. These solutions, given in section 1.3.2, are all 
basically unstable and shall not be discussed further. 



2.2 STABILITY OF THE PERIODIC SOLUTION 

The solution of equation 10 in section 1.3.2 yielded a 
periodic solution (equation 11) when e = -1. The solution 
shall be investigated for stability using the method of charac- 
teristic exponents. 

Suppose that one has a system such that 

du 



dt 



= P(u,v) , 



and 



dv 

dt 



Q(u,v) . 
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Suppose the solutions to these equations are known to be 



U 1 = u i^ t ^ and v i = 



Since these solutions are assumed to be periodic, they 
must have a period which shall be called T. 

Then define m 






• 3P 
•3u 






Then the periodic motion is stable if h > 0, and unstable 
if h < 0 [4]. Applying this approach to the specific problem 
at hand one obtains 



and 



d<j> 



dS 



= ^ 



dS 



= Q 

1 - 6 2 



These differential equations are essentially equivalent to 

derivatives with respect to time if the analysis is carried 

out for constant position, that is, 

d£ _ 1 d$ 

dS 3 dx 



at a given value of x since 

dS = dx - 3 dx . 

Carrying through the other operations one finds that 



and 



9P = = o 

3<f> 



3Q 3 

d\p dip x_g2 



o . 



Therefore h 0. 
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The theory then predicts that the motion is indifferent 
to perturbations. That is, the perturbation results in a new 
amplitude of oscillation which will be maintained. 

2.3 STABILITY OF THE PERTURBED STATIC FUNCTION 
2.3.1 THE WAVE EQUATION FOR THE PERTURBATION 

One can derive a wave equation for a perturbation applied 
to the static solution. As stated in section 1.2.3, only the 
case for e = -1 shall be analyzed so the static potential is 
represented by equation (7) . 

The wave equation of interest is 

V 2 <}> + A 2 4> 5 - - 2 = 0 . (12) 

c 9t 2 

Assume that 

4 > = 4> 0 + ^ , 

where \ p is very small. 

Since 




equation (12) becomes 

V 2 4> n + V 2 ip + A 2 (4> n +<M 5 - ^2 ^ = 0. 

0 0 C 3t 2 

The third term may be expanded to 

A 2 (4>q+54>q^+. • .+^ 5 ) 

But if i() is very small one may neglect all but the first two 
terms. Then, since 

V 2 <J> 0 + A 2 <J> q = 0 , 

one may write 

V 2 ip + 5A 2 cJ)^ - ^ = 0. (13) 

0 c at 2 
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Equation (13) shall be used as the wave equation for a 
perturbation applied to - 

This problem is of interest because, if <J>q reoresents any- 
thing like a particle, then the question arises whether the 
nonlinear self-interaction allows a stable particle or whether 
it tends to break up because of application of perturbations. 

2.3.2 LIAPOUNOFF'S STABILITY ANALYSIS AND STABILITY OF 
EQUILIBRIUM 

Suppose that one picks a point in space and proceeds to 
analyze the reaction of a perturbation approaching that point. 
Then <f>Q in equation (13) is a constant. Using the substitu- 
tions and the method outlined for the periodic solution in 
section 2.2, one obtains for a one dimensional perturbation 



is a constant. 

Since the equation was derived with x as a constant, the 
derivative is essentially with respect to time. One can now 
do a local stability analysis. Where the word local is used 
to emphasize the fact that only one point in space is being in- 
vestigated. The equation derived above obviously represents 
periodic motion if 3 2 < 1, <<nd an unbounded solution if 3 2 

However a nonlinear Wctve equation has been converted to 
a linear equation by assuming that ^ is small. One might ask 



d 2 \p 

dS 2 




(14) 



where 
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if such a prediction has any validity for the nonlinear problem. 
Liapounoff's stability analysis will yield the answer to this 
question . 

Let 



and 



du 

dt 



P(u,v) , 



and 



dv 

dt 



Q(u,v) , 



dv _ P _ cu+dv+Q 1 (u , v) 
du - Q au+bv+P 1 (u, v) ' 

where Q' and P' have no first order and no constant terms. If 

u and v are near zero, the last equation can be reduced to 

dv _ cu+dv 
du - au+bv ' 

and the equation has been linearized. 

Define s so that 



a-s c 
b d-s 



0 . 



(15) 



Liapounoff's stability criteria are as follows, if the 
roots Sj and S 2 of the determinant: 

(1) are complex or real, and Re(s^, S 2 ) are negative, 
then the solution is stable. 

(2) are complex or real, and at least one of the 
roots has a positive real part, the solution is unstable. 

(3) are imaginary and Re(s^, S 2 ) = 0, then the 
stability analysis is questionable. The question as to the 
solution's stability is not conclusively answered. 

Applying this analysis to equation (14) one proceeds 
as follows: 
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Let 



d^ _ - 
dS " ^ 



then 




and 



d£ _ aip 

HijJ 5 



Then in equation (15) 



c = a, b = 1, and d = a = 0 



Therefore s 2 -a = 0,s = ±/oi 



If 8 2 > 1, a > 0, and the roots of s are real 




Therefore, one can say that if 6 2 > 1, the solution, as 
already stated, is unstable. 

If B 2 <1, the roots of s are imaginary. Therefore, the 
analysis does not give reliable information concerning the 
stability of the motion. However, the motion is periodic and 
the method of characteristic exponents in section 2.2 mav be 
used. 

Applying this method once again 







= Q = 




Then 
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Therefore 



h 



T 



1 

T 



3P 

3 ^ 



+ 



3Q 

3 ? 



dt 



0 . 



Since h is zero, one concludes that the motion is indifferent 
to perturbations. This result does not change if we take into 
account all higher order terms in the expansion of ( q + ^) 5 , 
still considering 4 >q constant. The question arises at this 
point as to whether this approximation to "local" stability 
or instability, by keeping (|>q a constant and not really con- 
sidering the perturbation as propagating into regions of 
changing <j> Q , might invalidate the conclusions. It certainly 
is not giving information about propagating instabilities, 
known also as convective instability. In order to investigate 
this point the method in the next section will allow removal 
of the restriction that <J)q remain constant. The stability 
analysis technique to be used is well known in plasma physics. 



2.4 STABILITY ANALYSIS WITH THE WKB METHOD 

Another general approach to stability has been developed 
by Sturrock. It requires that the dispersion relation to the 
wave equation be obtained. For this analysis a three dimen- 
sional, spherically symmetric solution shall be obtained for a 
perturbation on the static potential function. 

2.4.1 AN APPROXIMATE SOLUTION TO THE LINEARIZED WAVE EQUATION 
In section 2.3.1, the linearized wave equation for a per- 
turbation on the static potential function was obtained: 

v 2 ip + 5X 2 <j>> - | 2 ill = 0 . 

0 C 3 1 2 
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If the static potential function (equation 7) is put in 
this equation, one obtains 

v 2 i|> + iSjsiS: — - i 2 111 = o . 

(r 2 +K 2 ) 2 c 3t 2 

It will be assumed that the perturbation , ip , is spher- 
ically symmetric, thus simplifying the mathematics involved, 
and amounts to considering s-states only. 

The equation then becomes 



3 2 i p + 2 _ 3jJj. + 15k 2 ^ _ 1_ 3 2 ip _ q 

3r 2 r 3r (r 2 +< 2 ) 2 c2 3 t 2 

To solve this equation we assume the solution to be of the 

form 



ip = T (t) R(r) , 



and a separation 

1 d 2 R + 
R dr 2 



constant a 2 . Then we mav write: 

2_ dR 15k 2 = _1 d^T = _ 2 

Rr dr (r 2 +K 2 ) 2 Tc2 dt 2 



d 2 T 

dt 2 



+ 



Ta 2 c 2 



0. 



T = e ±iut , » = ac . 



Sis + 1 ® + r \ 15>;2 + a 2 

dr 2 r dr l (r 2 « 2 ) 2 



= 0 



Let rR = u , then 

d 2 u 

dr 2 



+ u 



15 

(r 2 +< 2 ) 2 



+ a' 



= 0 



Let 



= kx, a 2 < 2 = Q 2 , 
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then 



d 2 u 
dx 2 



+ u 



15 



(x 2 +l) 



+ £ 2 ' 



0 . 



The exact solution to this equation is rather difficult 
to obtain. Therefore, the WKB approximation is used to get 
a solution which may be used for stability analysis. 

If one has a differential equation, 

+ K 2 u = o , 

dx 2 

then the WKB approximation is to use 



u 



K -i e i/Kdx 



For this case, 



K 2 



15 

(x 2 +l) 2 



+ ft 2 



r 



and 

K~% = (x 2 +l) % . 

[15+Q 2 (x 2 +l) 2 



The evaluation of the integral in the exponent will 
yield elliptic functions of all three tyoes . Such a func- 
tion will be very difficult to analyze for stability. There- 
fore, this integral will also be approximated. To do this 
the Q 2 , x plane shall be divided into two regions. One near 
the origin and one further from the origin. (See Figure 9) . 
(a) Assume that 

— — — < i , 

£2 2 (x 2 + 1 ) 2 

i.e., £2 2 (x 2 + l) 2 is large. 
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Then 



K = 



15 



(x 2 +l) 2 
K - fi 



4- fi ‘ 

r 



= n 



15 



Kdx - fix + 



1 + % 



15x 



Ifi 2 ( x 2 + 1 ) 2 

15 

fi 2 (x 2 +l) 2 



+ 1 



% 



4 fi (x 2 +l) 



15 . -1 

+ ?fi tan x . 



Then 



(b) Now assume that 



ft 2 ( x 2 + 1 ) 2 



15 



< 1 , 



i.e., ft 2 (x +1) is very small. 



K = 



/I~5 
x 2 + l 



1 + fi 



-1 



2 (X 2 +l) 

15 



T % 



Kdx =; /T5 tan ^x + 



( ix 3 + x ) . 

2/15 3 



To summarize: 
if 



15 



ft 2 (x 2 +l) 2 



< 1 , 



u = 



if 



u = 



(x 2 +l) % 



[15+fi 2 (x 2 + l) 2 



EXP i 



fix + 



15x , 15 , -1 

+ — tan x 



4fi (x 2 +l) 4fi 



fi 2 (x 2 +l) 2 



15 



< 1 , 



(x 2 + l) % 



[15+fi 2 (x 2 +l 2 ] % 



EXP i 



/15 tan ^x + 



fi‘ 



2/15 



(^x 3 +x) 



Quite a few approximations have been made in obtaining 
these solutions. Therefore, the regions in which the approx 
imations are useable must be determined. The condition for 
the validity of the WKB approximation shall be derived first 
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4 



ft 2 




Figure 9 
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In the differential equation. 



d 2 u 
dx 2 



+ K 2 u = 0 , 



make the independent variable substitution 

y = /Kdx . 

Then 



K 2 + K ' ~ + K 2 u = 0 . 
. 2 dy 



dy' 



where 



d 2 u K_^ du + u = 0 , 



dy' 



K 



2 dy 



K' = 



dK 

dx 



Now let 



u = vs , s = K ^e^ y 



ds ds 

dy ” K dx 



ds 



dy K 



= F -% 



K' 

3 

K 2 



+ iK' 



iy 



d 2 s 1 
,2 K 



dy ' 



5 (K ’ ) 



I \ 2 



7 - * — - 



K" iK' 



K 



Sh + 2ip- + v 

dy 2 d y 



K 

3 (K' ) 2 

4 K 4 



K 



- K' 



,!y 



— %• 



K 3 



= 0 . 



If the last term is zero, then the general solution for 



u is 

u _ K ~% (ae + be iy ) . 

Therefore the particular solution for u assumed by the WKB 
approximation is good when the absolute value of the last 
coefficient is very small. 
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If e is defined to be equal to this term, then by dif- 
ferentiating and by algebraic manipulation one obtains for 
this expression: 

_ 15 [15+fl 2 (x 2 +l) 2 ] - 75fl 2 x 2 (x 2 +l) 2 

[15 +Q 2 ( x 2 + 1 ) 2 ] 3 

The shaded area near the bottom of Figure 9 is the re- 
gion of space for which |e| > 0.1 . Within this region the 
WKB approximation is either not valid or of questionable 
accuracy . 

Next, consider the approximation of / Kdx. The approx- 

p 

imations used are those for the square root of (z+1) * . 

(z+1) * = 1 + %z - %z 2 ... , z < 1 



where only the first two terms were used in the expansion. 

The error is not larger than the absolute value of the 
third term. Again, the value, 0.1, is chosen. For case a, 

15 2 



k -z. = 



< .1 



4Q 2 (x 2 + 1 ) 



Making the expression an equality, we can write 
ft 2 (x 2 +l) 2 = 23.75 , or roughly , 

n 2 = -2! — . 

(x 2 +l) 2 

If the same procedure is followed for case b, one obtains 

ft 2 = . 

(x 2 + 1 ) 2 

These two functions are plotted in Figure 9 . The shaded 
region between them represents the region where the approx- 
imation to the integral is unreliable. The region near the 
origin is represented by case b; the region above the shaded 
area by case a. 
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2 . 4 „ 2 STURROCK’S STABILITY ANALYSIS 



Suppose that the solution of a wave equation yields 

, i(kx-iot) 

<P = e , 

and 

k = k(cj) , 
to = to(k) . 

If to is plotted in the uk , to , plane for all real k, then 
the shape of the plot will yield information about the stabil- 
ity of the wave. Furthermore, a plot of k in the k^, k^ plane 
for all real to will yield information concerning the type of 
instability of the wave. 

For purposes of clarity the possible alternatives shall 
be divided into four cases. 

Case (a) -evanescent waves 

These waves are stable. The system is stable if 
for real k, to is always real. This situation is diagrammed 
in Figure 10. 





0) . 




l 


' \ 


/ 0) 

C m — " - 


— — <■ > 





FIGURE 10 

Dispersion Relation in the Complex to Plane 
for Evanescent Wave 
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Case (b) - non-convecti ve instability 

This case represents waves in which a disturbance 
grows and expands about the point at which it occurs , but does 
not travel. 

Such a case occurs in the type of plots diagrammed in 
Figures 11 and 12. Note that co^/w r -> 0 as k -*■ ± <*> , and that 
k is real for all real w. 




Dispersion Relation in the Complex co Plane 
for a Non-convecti ve Instability 




FIGURE 12 

Dispersion Relation in the Complex k Plane 
for a Non-convective Instability 
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Case (c) - convective instability 

This case represents a wave disturbance which qrows 
and propagates. The ur , u> r plot is the same as Figure 11, 
but the plot now also has complex values of k for real to and 
k./k 0 as to i ». The plot of k. , k would be the same 
as Figure 11 with the appropriate changes in labeling of 
axes . 

Case (d) In this case ur/io^. -*■ constant ^ 0 as k ± 00 . 
This condition is represented in Figure 13. 




Dispersion Relation in the Complex to Plane 
for an Undefined Instability 



The only conclusion that can be drawn in this case is 
that the disturbance is unstable. It is Sturrock's conten- 
tion that this situation will occur only when the model is 
inexact or incomplete. 

To apply Sturrock's theory to the WKB solution obtained 
for the perturbed wave function, let 



k = ft + X - ?— + -H tan 1 x + i- In — 

4£2 (x 2 + l ) 4Qx (x 2 + l) 



+ Q ‘ 



(16) 
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when 



15 



< 1 ; 



ft 2 ( x 2 +1 ) 2 



and 



when 



k = tan ^ x + — ( ix 2 + 1 ) + 4— In 

x 2/15 3 4x 

ft 2 (x 2 + l) 2 



+ ft 2 



(x 2 +l) 2 



(17) 



15 



< 1 . 



Plots of ft. vs. ft and k. vs. k are shown in Figures 14 
l r i r ’ 

and 15 for x =1.5. The plot of k^ , k^ was obtained directly 
from equations (16) and (17) . The plot of ft^, ft^ was obtained 
by approximating the logarithms and then solving for ft . 

For equation (16) at x = 1 . 5 , 

k = ft + iili + i ( . 166) In [1. 42+ft 2 ] . (18) 

The logarithm is approximated by factoring out ft 2 and 
dropping the remainder, 

In [1 . 42+ft 2 ] = 21n ft + In + 1] . 

ft 2 

But from the restrictions on equation (18) , due to approxima- 
tions in the last section, the last term cannot be greater 
than In 1.63. Since for this equation |ft| >.1.51, the last 
term may be dropped if ft is very large. The logarithm of ft 
is then approximated by 



Although this introduces a great deal of error if ft be- 
comes too large, it does allow one to solve for ft using only 
a quadratic and the curve is correct in general features. 

The result is 

k = ft + — ^ — + i(.33 3) — , 
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Dispersion Relation in Complex oo Plane for the 
Perturbation of the Static Solution 



Figure 14 




Dispersion Relation in Complex k Plane for the 
Perturbation of the Static Solution 

Figure 15 
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which may be solved for ft, 



ft = | - i ( . 166 ) ± %"\/k 2 -14.56+|i (2-k) . (19) 

The same general ideas may be carried out for equation 

(17) . 

k = 2.54 + .266ft 2 + i ( . 117) ft 2 + i(.0585) 

ft - ± 2.04e -1 11,90 ° Vk-2.54-i(.0585) (20) 

This equation is useful only if |ft| < .95. This restric- 
tion severely limits the values allowed by k and the equation 
is not very useful in obtaining Figure 14. The plot was ob- 
tained using equation (19). When the ft^ , ft r plot is dis- 
placed below the ft r axis, the motion is stable. This is the 
case for large values of ft r in Figure 14. Nearer the origin 
one can see that the plot looks as though it is unstable. In 
this region ^ > 1 . This corresponds to 6 > 1 in section 
2.3.2 where an unstable solution was predicted. The upper 
curve in this region has been terminated because neither equa- 
tion 19 nor 20 yield valid values of ft. However, the curve 
seems to curve back towards the origin from calculations with 
equation (19) . This would be consistent with the idea that 
the solution becomes stable for large values of ft. 

The k^, plot, Figure 15, does not give as much infor- 
mation. The large gap between the two line segments is due 
to the limitations on equations (19) and (20) . 

It can be shown with equation (16) that 

. k . 

lxm l _ /. 

ft-* 00 k 

r 
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However,, due to the limitations on the ranges of equations 
(16) and (17) , one cannot determine very much about the aaos 
in Figures 14 and 15. Therefore, one cannot say much about 
the type of instability that exists for small k. 
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CHAPTER III 



SUMMARY AND CONCLUSIONS 

3 . 1 SUMMARY 

The nonlinear wave equation, 

V 2 4> - eX 2 <j) 5 - i-2 ** 1 = 0 , 

c at 2 

derived from the Klem-Gordon wave equation has been investi- 
gated. The e = +1, this equation represents a true analog 
to the Klein-Gordon wave equation with the mass replaced by 
field self-interaction. For e = -1, the nonlinear term may 
be interpreted as the source term in the electromagnetic 
potential equation with the source given by field self- 
interaction . 

Solutions to the static wave equation are characterized 
by a fundamental length k and an "energy parameter", h, as 
integration constants. The term "energy parameter" is used 
only in the mechanical analog sense. It is found that for 
£ = +1, the solutions all are unbounded and not everywhere 
real for all values of h. The case h = 0 may behave either 
as a 1/r potential, infinite at the fundamental distance k, 
and with a nonreal solution for r < K ; or as a potential well, 
centered at the origin, and infinite at r = k with no real 
solution for r > k. The first result may allow the inter- 
pretation of a particle with a "hard core". The second re- 
sult could be interpreted as the basis of a particle whose 
characteristics are determined by eigenfunctions within the 
potential well. 
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For e = -1,, the solutions are everywhere real. In partic- 
ular, the solution with h = 0 is also everywhere finite and 
has a 1/r dependence for r > k. This solution may represent 
an electromagnetic model of a charge with k as the extension 
of the charge. 

Solutions of the full wave equation were obtained. The 
solutions are distinguished in their character by their veloc- 
ity relative to light and by the sign of e. Those solutions 
for which 8 < 1, e = 1 or M 1, e = -1, are unbounded and 
not everywhere real. Those solutions for which 8 > 1, e = 1, 
or 8 < 1, e = -1 are periodic. 

The study of the stability of the solutions with £ = -1 
represented the major objective of this paper. The stability 
was tested in three different ways. The stability of the 
periodic solutions was tested for one dimensional solutions 
using the method of characteristic exponents used in non- 
linear mechanics. These solutions were found to be indif- 
ferent to perturbations, that is the perturbation results in 
a new solution which will be maintained without blowing up or 
returning to the original. 

Next, the stability of a perturbation on the static po- 
tential for h = 0 was investigated. The wave equation for 
the perturbation was linearized and analyzed for stability 
assuming that <j> Q is constant. This amounts to studying "local 
stability". If B 2 <1, Liapounoff's stability analysis gives 
no information about the stability of equilibrium. The per- 
turbation xn this case is periodic and may be analyzed for 
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stability of the focal point of the phase diagram using the 
method of characteristic exponents. The result found was 
again indifferent to stability. If B 2 > 1, Liapounoff's 
method indicates instability. 

Finally, the linearized perturbation equation was inves- 
tigated using Sturrock's analysis of the dispersion relation 
obtained from a WKB approximation to the solution. The dis- 
persion relation obtained was investigated at r/< = 1.5, 
i.e., near the center of the charge. Of the two possible 
complex roots of w for real wave number, k, one was found to 
have only negative imaginary parts indicating a stable solu- 
tion. The other solution has positive imaginary values 
around w r = 0 , indicating instability when $ > 1. The char- 
acter of the instability could not be determined because of 
restrictions on the equations used due to the WKB and other 
approximations . 

The results obtained from Sturrock's analysis agrees with 
the results obtained from the first two methods, if 8 > 1. 

For 8 < 1, Sturrock's method predicts stability, whereas the 
method of characteristic exponents predicts indifference. 

3.2 CONCLUSIONS 

It has been demonstrated in the paper that a relatively 
simple nonlinear field equation can serve as a convenient test 
ground for investigation of stability concerning solutions of 
nonlinear field equations which might represent particles or 
electromagnetic sources arising from field self- interaction . 
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The analysis in this case has shown that certain perturba- 
tions are indifferent or stable. The result seems to be sen- 
sitive to the change from consideration of "local stability" 
to consideration of the varying <J>q. This suggests that a non- 
linear stability analysis might eventually change the result 
even further. To determine whether or not this is true will 
require a better approximation to the solution for the per- 
turbation. This may be accomplished by including nonlinear 
terms in the linearized equation and by carry inq the WKB 
method to high orders of approximation by some kind of itera- 
tion. A second alternative would be to seek new avenues of 
stability analysis, eventually using computer techniques. 

It has been shown that there exist other classes of per- 
turbations which are unbounded. It remains to be established 
whether or not these perturbations are everywhere real and 
what their physical interpretation would be, perhaps connected 
with shock-like phenomenon. Answers to these questions again 
lie in improvement of the WKB approximation and the analysis 
of the dispersion relation. 
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APPENDIX A 



The General Solution of the Static Wave Equation 

In the main body of the paper potential functions were 
obtained for the case h = 0 in the equation 

jdFj 2 = 2h + £ ^ 2f6 + %F 2 . (4) 

In general, the solutions of the equation are elliptic 
functions whose form is dependent upon the sign of h and e, 
and upon the relative size of h and A. 

In equation (4) , let F = v 2 and 

24h = Sa 2 

where 

6 = ±1 . 



The equation then becomes 



3 ( = vU4A 2 v 3 + 3v+<5a 2 ) 



This equation, when it yields elliptic functions, may be 
evaluated using a handbook of elliptic functions. [6] The 
cubic on the right must first be factored. 



4 A ' 



dv 

du 



r , , n \ i r . A+B (A— B) , — — i r A+B (A— B) , — =j-i ^ 

e v{ v- (A+B) ] [v + — ~ — * Z-3] [v + ~ /-3J } 



( 21 ) 



where 



A = 



B = 



2 A 



2 A 



<5a^A 



6a 2 A 



+ V 



a 1 * A 2 + e 
o 



- VajA 2 +e 



( 22 ) 



(23) 
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To solve the equation, the solution shall be divided into two 
cases . 

Case 1. e = 1 

If 6 = 1 (i.e., h > 0) then F can have all values and 

0 £ v < «> . 

The differential equation is 




dv 



0 [v ( v-R 1 ) (v-R 2 ) (v-R 3 ) ] * 






where R has been substituted for the roots, 
A and B are 



For this case, 



51 t- a 0 X Wa^+l] 3 



From an elliptic functions handbook the solution to this equa- 
tion is obtained by letting 

2 _ (a+8) 1 - (a+b) 2 



k = 



4a8 



Q _ 2X (a-8 ) v- (a+b) a 
COS " 2X (a+8) v- (a+b) a ' 



Then 



a = ( a 2 + 1 + b 2 ) % 

8 = [ 3 (a 2 - 1 + b 2 ) ] % 
a = 2XA, b = 2XB 



Q fa8 n r 

cos 0 = cn In — , 

3 K 




where the modulus of the elliptic function is k. 
Let the argument of the cn be x. Then 



*0 = 



| v \h _ j i 


a (a+b) tn 2 x/2 dn 2 x/2 


lr/ ’ 2Xr 


,atn 2 x/2 dn 2 x/2-8 



A 



% 
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This solution actually gives a series of unbounded waves. 
The complicated function is real for small values of 

tn 2 x/2 dn 2 x/2 

because a + b is negative. 

If 6 = -1 (i.e., h < 0 the differential equation is 



__ dv , _ 1r , r 

Tv(v-R ] TTv-R 2 ) (v-R 3 ) P ~ ln k 

R 1 

where the roots are defined by equations (21) , (22) , and 

(23). Specifically A and B are 

1 

a = 2AA = [a 2 A +V a^X 2 +l] 3 

b = 2AB = [a 2 A - Va*A 2 +l] 

and R^ = A + B 

The solution will be very similar to that for 6=1, 
except that from section 1.2.1, F has a minimum of 

F min = (A+B ) ^ 
or v min = A + B. 

This will necessitate redefining slightly the elliptic func- 
tion and its modulus. 

i, 2 _ (a+B) 2 - (a+b) 2 
K 

a 2 = 3(a 2 -l+b 2 ) 



/3 
2 A 



B 2 = a 2 + 1 + b 2 



cos 0 



2 A (a-B ) v + (a+b) B 

2 A (a+B ) v - (a+b) B 
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„ / a8 , r 

cos 0 = cn , — In — = cn x 



f 1 


6 (a+b) \ 


1 2 Ar 


. 8 -atn 2 x /2 dn 2 x/ 2 j 



This function also represents an unbounded wave. The wave 
minimum occurs periodically and approaches zero as r increases. 



Case 2. e = -1 

If 6 = 1, the differential equation will be 



4 A : 



t^j = -v[(v-R 1 ) (v-R 2 ) (v-R 3 ) ] 



whe 



re R^ , R 2 » R 3 are defined by equations ( 21 ) , ( 22 ) , and 



(23) . A and B will be 



a = 2AA = [a 2 A + v^gA. 2 -!] 



b = 2 AB = [a 2 A -VajA 2 -!] 3 . 



In this instance there are three different solutions, 
that 



a o A ' 



> 1 



Then 



/3 
2 A 



- v 



dv 



[ v(R^-v) ( v-R 2 ) (v-R 3 ) ] 2 



k = 



In - 

K 



where the roots will be specified as 
R^ = A + B and 

R 2 and R 3 are the complex roots. 
From section 1.2.2, V has the range 

0 < v < A + B . 



Assume 
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Let 



k 



2 _ (a+b) 2 - (a-B) 2 

2aB 




a 



2 



3 (a 2 +l+b 2 ) , B 2 = a 2 - 1 + b 2 



Then 



cos 6 = cn 




cn x , 



where the elliptic function has the modulus k. Prom the last 
equation one may solve for v and <t> Q , 



<t>Q oscillates in this case with decreasing amplitude as r in- 
creases . 

Assume a^A 2 = 1 • This is a special case in which the 
solution may be expressed in terms of elementary functions. 
One has 



1 / B (a+b) tn 2 x/2dn 2 x/2 



,^ r a+Btn 2 x/2dn 2 x/2 



/3 

2 A 




dv 




J 0 



Make the substitution 



1 




Then from a table of integrals : 




y 



/3 
2 A 



h\ 



Upon transposing the sin 1 and some manipulation, 

- - I j-~ cos ln r/< | ^ 

^0 _ 2Ar 1 2+cos In r/K I 
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Assume 



a^X 2 >1. A better representation of the roots 
of the cubic in this case would be 

^ cos <j>/3 , 

R 2 = ~ 2 j (cos 4>/3+/3sin 4>/3 ) , 

R 3 = (-cos 4>/ 3+/3sin <p/3) , 

cos (p = a^X . 

The range of v as calculated in section 1.2.2 is 
0 < v < T- cos 4>/3 . 



The differential equation is 

-V 

/ 3 ' 

2 X 



dv 



k 



- In 



1/1 cos V3 ( V -R 2 > (v-R 3 > ] " 



2 sin 



2<j) 



Define : 



k 2 = 



ns 2 (p ^ 2 <p 

/3 cos — j + sin — j 



Then 



Q _ ( [cos <j>/3+/Tsin 4>/3 ] [cos <j)/3-Xv] 
s 1 n w _ v 

\ cos (p/3 [ 2 Xv+cos <p/ 3+/3sin <p/ 3] 



Q ,1/ 3cos 2cf>/3+sin 2<j)/3]^ 

sin 9 = sn % — — — 

1 /3 / 



i K 
In — 
r 



or 



sin 0 = sn x 



From the last equation v and q may be obtained, 



*0 = 



_1 

Xr 



cn 2 Xcos <j>/3(cos cj)/3+/3sin 4>/3) 
.cos (j>/3 (l+2sn ? X ) + /3sin (p/3 



h 



Let 6 = -1. Then the differential equation is 

_ 2 



3 rdvr 
4X 2 [duj " v 



3v 



0 



v “ 4X 2 + IT 2 
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7S|hi 



A and B in the roots as defined by equations (21) , (22) , and 

(23) are ^ 

a = 2XA = [-a 2 X + \ajx 2 -l] 3 

1 

b = 2XB = [-a 2 X - VaJX 2 -l]^ 

It was shown in section 1.2.2 that when h < 0, F cannot 
be real unless, a^X 2 < 1“ Therefore this will be the only 
case considered. The roots can be better expressed as 

= j cos 4/3 , 

R 2 = 2X ^ -cos 4>/3-/3sin 4>/3 ) , 

= 2 y (-cos 4/3+v / 3sin 4>/3 ) , 

cos 4> = -a 2 X , 



and v has the range 



Rg v _< R^ 



Therefore 



/3 

2X 



- v 



R, 



dv 

[v (R, -v) ( v-R~ ) (v-R,) ] 



k = In — 



Define 



2 _ /3cos 24/3+sin 24/3 
~ 2sin 24/3 



sin 6 = / 4/3sin 4/3 (co s 4/3-Xv) 

V (3cos 4/3-»^3sin 4/3) (2Xv+cos 4/3+/3sin 4/3) 



Then 



sin 0 = sn 



sin 



„24/3j 



h 



In — 
r 



= sn x 
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where the sn x has modulus k. 
solving for v one may obtain 



By squaring both sides and 



*0 = 



ri 


sin 2<f>/3 (2-sn" x) -/3sn 2 xcos 


24>/3\l 


r 


sin <f>/3 (2-sin 2 x) +/3sn 2 xcos 


4>/3 }j 



h 



This is an oscillating function. In the limit as cos 9 
the term in the brackets becomes identially 1. Then 



*0 = 



2 Ar 



% 



- 1 , 



and F is a constant as predicted in section 1.2.2. 
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APPENDIX B 



Exact Solution of /Kdx 

The exact solution of this integral involves all three 
types of elliptic functions. One may obtain the solution in 

the following manner. 

Kdx = U5+ft 2 (x 2 + l) 2 ] % dx 
x 2 + l 

Let z = x 2 + 1. Then the integral becomes 

[15+ft 2 z 2 ] ^dz 
2z(z-l)* 

This may be reduced to 



15dz 


ft z dz 


2ftz(z-l)^ (— + z 2 \^ 


2(z-l) % + z 2 | % 


U 2 J 


U 2 1 



De fine : 



0 = (15+ft 2 )^ , 

2 _ 8 - £2 
" 20 



cos 0 



0 + ft-ftz 
0 - Si+fiz 



$ - fix 2 
0 + ftx 2 



1 



v = F(e,k) , 

where F(0,k) is the elliptic integral of the first kind. 
The integrand may be replaced by 

15 ^ (1+cn v)dv dv ^ [ B+ft+ (fl-8 ) cn v ] 

2 0 $+ft+(ft-B)cn v 2 0' 1+cn v 



where cn v is an elliptic function of v. 
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Integrating over v one obtains 

v|§|* - (B-a)*E(v) + ( 6 B)* sn vdn v 



1+cn v 



15 ie+n\ 
13-nl 



4 m) 



, 0 -3k 4 

T^O, , k) 



+ - — tan 



-1 



15 
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2 _ 

3ft 



sd v 



where E(v) is the elliptic integral of the second kind, 



TT (0 , 



-6k; 

n 



,k) is the elliptic integral of the third kind, and 



the functions sn v, dn v, sd v, are elliptic functions of v. 

The following equalities apply: 
sn v = sin 0 , 
cn v = cos 0 , 
dn v = (l-k 2 sin 2 0 ) % , 



sd v = 



sin 0 



(l-k 2 sin 2 0) 5 

Using these equalities the elliptic functions may be replaced 
by functions of x. The results yield: 



v (§} Mtt “ (3ft)^E(v) - 



B-n 



15 

4(ft3) % 



3+ft 



/fl 3k 4 . v 
it (0,-— ,k) 



+ -£ tan" 1 

/2 



15 x 



[B 2 +2ft 2 x 2 +ft 2 x 4 ]* 



ft 



B+ftx‘ 



(B 2 +2ft 2 x 2 +ft 2 x 4 ) 



k 
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